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ABSTRACT 

We investigate the lepton flavor violation in the framework of the MSSM with 
right-handed neutrinos taking the large mixing angle MSW solution in the quasi¬ 
degenerate and the inverse-hierarchical neutrino masses. We predict the branching 
ratio of fi —> e + 7 and r —> p + 7 processes assuming the degenerate right-handed 
Majorana neutrino masses. We find that the branching ratio in the quasi-degenerate 
neutrino mass spectrum is 100 times smaller than the ones in the inverse-hierarchical 
and the hierarchical neutrino spectra. We emphasize that the magnitude of t/ e3 is one 
of important ingredients to predict BR(p —> e + 7). The effect of the deviation from 
the complete-degenerate right-handed Majorana neutrino masses are also estimated. 
Furtheremore, we examine the S$ L x S^ R model, which gives the quasi-degenerate neu¬ 
trino masses, and the Shafi-Tavartkiladze model, which gives the inverse-hierarchical 
neutrino masses. Both predicted branching ratios of p —> e + 7 are smaller than the 
experimantal bound. 
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1 Introduction 


Super-Kamiokande has almost confirmed the neutrino oscillation in the atmospheric 
neutrinos, which favors the ry —> v T process |l]|. For the solar neutrinos [^J, [|, the 
recent data of the Super-Kamiokande and the SNO also suggest the neutrino oscillation 
v e —> z/ x with the large mixing angle (LMA) MSW solution, although other solutions 
are still allowed [Q. |S[]. If we take the LMA-MSW solution, neutrinos are massive and 
the flavor mixings are almost bi-maximal in the lepton sector. 

If neutrinos are massive and mixed in the SM, there exists a source of the lepton fla¬ 
vor violation (LFV) through the off-diagonal elements of the neutrino Yukawa coupling 
matrix. However, due to the smallness of the neutrino masses, the predicted branching 
ratios for these processes are so tiny that they are completely unobservable f 6 |. 

On the other hand, in the supersymmetric framework the situation is quite different. 
Many authors have already studied the LFV in the minimal supersymmetric standard 
model (MSSM) with right-handed neutrinos assuming the relevant neutrino mass ma¬ 
trix 0 . g. 0 . 0 . In the MSSM with soft breaking terms, there exist lepton flavor 
violating terms such as off-diagonal elements of slepton mass matrices (mf R )ij 

and trilinear couplings Af^-. Strong bounds on these matrix elements come from requir¬ 
ing branching ratios for LFV processes to be below observed ratios. For the present, the 
most stringent bound comes from the n —> e + 7 decay (BR(/i —> e + 7 ) < 1.2 x ICY 11 ) 


11], However, if the LFV occurs at tree level in the soft breaking terms, the branch¬ 


ing ratio of this process exceeds the experimental bound considerably. Therefore one 
assumes that the LFV does not occur at tree level in the soft parameters. This is 
realized by taking the assumption that soft parameters such as (m|)^., (m l R )ij, A®-, 
are universal i.e., proportional to the unit matrix. This assumption follows from the 
minimal supergravity (m-SUGRA). However, even though there is no flavor violation at 
tree level, it is generated by the effect of the renormalization group equations (RGE’s) 
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via neutrino Yukawa couplings. Suppose that neutrino masses are produced by the 
see-saw mechanism a, there are the right-handed neutrinos above a scale M R . Then 
neutrinos have the Yukawa coupling matrix Y v with off-diagonal entries in the basis of 
the diagonal charged-lepton Yukawa couplings. The off-diagonal elements of Y u drive 
off-diagonal ones in the and A- matrices through the RGE’s running | 13 [ . 

One can construct Y v by the recent data of neutrino oscillations. Assuming that 
oscillations need only accounting for the solar and the atmospheric neutrino data, we 
take the LMA-MSW solution of the solar neutrino. Then, the lepton mixing matrix, 
which may be called the MNS matrix or the MNSP matrix is given in ref. (T^j . 

Since the data of neutrino oscillations only indicate the differences of the mass square 
A mlp neutrinos have three possible mass spectra: the hierarchical spectrum m u3 
m v 2 m u 1 , the quasi-degenerate one m v i ~ m V 2 — m u3 and the inverse-hierarchical 
one m v i ~ rn u2 > m v3 . 

We have already analyzed the effect of neutrino Yukawa couplings for the /1 —> e + 7 


process assuming the quasi-degenerate one and the inverse-hierarchical one [|I 7 |. In this 
paper, we present detailed formula in our calculations of /i —> e + 7 and discuss the 
dependence of the SUSY breaking parameters for the branching ratio. In the previous 
paper, the right-handed Majorana neutrino masses are assumed to be completely de¬ 
generate. We study the effect of the deviation from this degeneracy in this work. The 
correlation between BR(/r —> e + 7) and BR(r —> /i + 7) is also calculated. Furthere- 
more, two specific models of the neutrino mass matrix are examined in the /j, —> e + 7 
process. 

This paper is organized as follows. In section 2 , we give the general form of Y v and 
YJ Y u , which play a crucial role in generating the LFV through the RGE’s running. In 
section 3 , we calculate the branching ratio of the processes /r —> e + 7 and r —> /1 + 7, 
respectively in the three neutrino mass spectra. In section 4 , we examine the S 3l x S 3r 
model, which gives the quasi-degenerate neutrino masses, and the Shafi-Tavartkiladze 
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model, which gives the inverse-hierarchical neutrino masses. In section 5, we summarize 
our results and give discussions. 


2 LFV in the MSSM with Right-handed Neutrinos 

2.1 Yukawa Couplings 


In this section, we introduce the general expression of the neutrino Yukawa coupling 
Y„, which is useful in the following arguments, and investigate the LFV triggered by 
the neutrino Yukawa couplings. The superpotential of the lepton sector is described as 
follows: 

MW = Y e LH d e c R + Y„LH u v c R + ^v c R T M R v c R , (2.1) 

where H U1 H d are chiral superfields for Higgs doublets, L is the left-handed lepton 
doublet, e R and v R are the right-handed charged lepton and the neutrino superfields, 
respectively. The Y e is the Yukawa coupling matrix for the charged lepton, M R is 
Majorana mass matrix of the right-handed neutrinos. We take Y e and M R to be 
diagonal. 

It is well-known that the neutrino mass matrix is given as 

m„ = (Y u v u ) T M r 1 (Y„u u ) , (2.2) 


via the see-saw mechanism, where v u is the vacuum expectation value (VEV) of Higgs 
H u . In eq.( [2.2|) , the Majorana mass term for left-handed neutrinos is not included since 
we consider the minimal extension of the MSSM. 

The neutrino mass matrix m„ is diagonalized by a single unitary matrix 


m 


diag _ TjT 


= U., NIO m,,U 


MNS 1AA F ^ MNS ? 


(2.3) 


where U MNS is the lepton mixing matrix. Following the expression in ref. [[Uj , we write 
the neutrino Yukawa coupling as 


Y v = — vW 8 R J mf; ag U 


MNS ’ 


(2.4) 
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or explicitly 


1 

(y/M^ 

0 

0 > 


( y/rfhX 

0 

0 \ T 


Y„ = - 

0 

V m r2 

0 

R 

0 

y/m v 2 

0 u MNS , 

(2.6) 

Vu 

l 0 

0 

V j 


l 0 

0 

y / rn^J 



where R is a 3 x 3 orthogonal matrix, which depends on models. Details are given in 
Appendix A. 

At first, let us take the degenerate right-handed Majorana masses M n \ = M r2 = 
M r 3 = Mr. This assumption is reasonable for the case of the quasi-degenerate neu¬ 
trino masses. Otherwise a big conspiracy would be needed between Y u and M R . This 
assumption is also taken for cases of the inverse-hierarchical and the hierarchical neu¬ 
trino masses. We also discuss later the effect of the deviation from the degenerate 
right-handed Majonara neutrino masses. 

Then, we get 


and 




Vm~ r r 


( 


0 

V o 


0 

y/m u 2 

0 


0 ^ 

0 U^s 
y/m^J 


Yj,Y„ 


or equivalently, 




U 


MNS 


( m u i 


0 

V o 


0 0 \ 

o u; NS , 

0 m„ 3 j 


( 2 . 6 ) 


(2.7) 


M ^ 

( Y i Y -L = . 


( 2 . 8 ) 


"7t i =1 


where U a p s are the elements of U MN s- It is remarked that Y^Y^ is independent of R 
in the case of M R \ = M r2 = M R 3 = M R . ft may be important to consider the deviation 
from the degenerate right-handed Majonara neutrino masses. Detailed discussions are 
given in subsection 3.2. 

Note that this representation of the Yukawa coupling is given at the electroweak 


scale. Since we need the Yukawa coupling at the GUT scale, eq. (|2.5|) should be modified 
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by taking account of the effect of the RGE’s 
at a scale M R are given as 


18, 19, 20]. Modified Yukawa couplings 


\fMZ 

(y/rrh ,1 

0 

0 ^ 


(l 

0 

0 ^ 


Y u = — —- R 

0 

sjm v 2 

0 

tJ M Ns y 

0 

1 

0 ’ 

(2.9) 

Vu 

l 0 

0 

yjm v 3 ) 


lo 

0 

Vh) 



with 



1 ftR 




1 ftR „ I 

I g = exp 

CiS,dt . 

, I t = exp 

v ‘ dt ] 

, I T = exp 

L H 


( 2 . 10 ) 


where t R = In M n and t z = lnM z . Here, gis(i = 1,2) are gauge couplings and y t 
and y r are Yukawa couplings, Cj’s are the constants (|, 3). We shall calculate the LFV 
numerically by using the modified Yukawa coupling in the following sections. 

As mentioned in the previous section, there are three possible neutrino mass spectra. 
The hierarchical type (m u i <C m v2 m„ 3 ) gives the neutrino mass spectrum as 


m u i ~ 0 , 


m u 2 = 



m y 3 = \j Ami 


atm ? 


( 2 . 11 ) 


the quasi-degenerate type {m v \ ~ m V 2 ~ m u3 ) gives 


1 

m v i = m v , m v2 = m v + -— 

2m u 


and the inverse-hierarchical type (m u i 


Am| , m v3 = m u + 
~ m v2 » m v3 ) gives 


1 

2 m u 


Am 


2 

atm ? 


( 2 . 12 ) 


= J AmL , mj,! = m u 2 ~ - -, m v3 ~ 0 . (2.13) 

v 2m u2 

We take the typical values Am 2 tm = 3 x 10 _ 3 eV 2 and Amg = 7 x 10 _ 5 eV 2 in our 
calculation of the LFV. 


We take the typical mixing angles of the LMA-MSW solution such as S 23 = l/\/2 


and s 12 = 0.6 [|16J, in which the lepton mixing matrix is given in terms of the standard 


parametrization of the mixing matrix |23 


as follows: 


U 


MNS 


1 G3C12 C13S12 5 i 3 e"^\ 

-C 23 S 12 - s 23 s 13 c 12 e 1 ' 1 ’ C 23 C 12 - s 23 s 13 si 2 e lrt> S 23 C 13 , 

V S 23 S 12 - c 23 si 3 Ci 2 e^ -S 23 C 12 - c 23 s 13 s 12 e l(1, c 23 ci 3 ) 


(2.14) 
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where s t] = sin d tJ and Cij = cos 9{j are mixings in vacuum, and 0 is the CP violating 
phase. The reacter experiment of CHOOZ fU] presented a upper bound of s 13 . We use 
the constraint from the two flavor analysis, which is S 13 < 0.2 in our calculation. If we 
take account of the recent result of the three flavor analysis | 22 ], the upper bound of 
S 13 may be smaller than 0.2. Then, if we use the results in |22[], our results of /1 —* e + 7 
are reduced at most by a factor of two. In our calculation, the CP violating phase is 
neglected for simplicity. 


2.2 LFV in Slepton Masses 

Since SUSY is spontaneously broken at the low energy, we consider the MSSM with 
the soft SUSY breaking terms: 

- £ soft = (m DijQlQj + (m l)ijU* Ri u Rj + 

+ ( m z) ijULj + (m DijeCe-Rj + (m l)ijC Ri D Rj 

+™H d H d H d + m 2 Hu HIH u + (BfiH d H u + ^B uij M Rij v* Ri p* Rj + h.c .) 

+(A d ij H d d* Rj q j + A u ij H u n* Ri Q j + + A 

+ l -M x BlBl + hd 2 W£W a L + hhG a G a + h.c.) , ( 2 . 15 ) 

where mj~, m|, m^, m|, m? and m? are mass-squares of the left-handed squark, the 
right-handed up squark, the right-handed down squark, the left-handed charged slep¬ 
ton, the right-handed charged slepton and the sneutrino, respectively. The and 
rh 2 u u are mass-squares of Higgs, A d , A u , A e and A u are A-parameters for squarks and 
sleptons, and Mi,M 2 and M 3 are the gaugino masses, respectively. 

Note that the lepton flavor violating processes come from diagrams including non¬ 
zero off-diagonal elements of the soft parameter. I 11 this paper we assume the m- 
SUGRA, therefore we put the assumption of universality for soft SUSY breaking terms 
at the unification scale: 

( m zk = ( m !)u = = ■■■ = Sijml , 
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m H d = m H u = m 0 , 

A u = Y u a 0 m 0l A e = Y e a 0 m 0 , 

A u = Y u a 0 m 0 , A d = Y d a 0 m 0 , (2-16) 

where mo and a 0 stand for the universal scalar mass and the universal A-parameter, 
respectively. Because of the universality, the LFV is not caused at the unification scale. 

To estimate the soft parameters at the low energy, we need to know the effect of 
radiative corrections. As a result, the lepton flavor conservation is violated at the low 
energy. 

The RGE’s for the left-handed slcpton soft mass are given by 
= h-7-( m zk 

“P MSSM 

+ Ieb + Y I Y - m i)« + 2 < Y I m ^ Y - + *k Y IW + a'a„)«] , 

(2.17) 

while the hrst term in the right hand side is the normal MSSM term which has no 
LFV, and the second one is a source of the LFV through the off-diagonal elements of 
neutrino Yukawa couplings. The RGE’s are summarized in Appendix B. 

3 Numerical Analyses of Branching Ratios 

Let us calculate the branching ratio of e* —> ej + 7 (j < i). The amplitude of this 
process is given as 

T = ee a *(q)u j (p)m ei ia a pq l \A L P L + A R P R )ui(q - p ) , (3.1) 

where Ui is the wave function of i-th charged lepton e*, p and q are momenta of ej 
and photon, respectively, e is the electric charge, e is the polarization vector of photon, 
and P LJt are projection operators : P L , R = (1 =F 7 s)/ 2 . The A L,R are decay amplitudes 
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Figure 1: Feynman diagrams which contribute to the branching ratio of e* —> ej + 7 . 
There are two types of diagrams, (a) neutralino-slepton loop and (b) chargino-sneutrino 
loop. 


and explicit forms are given in Appendix C. It is easy to see that this process changes 
chirality of the charged lepton. The decay rate can be calculated using A L,R as 

r(e, - ej + 7 ) = -^<(\ AL \ 2 + I^T) • (3.2) 

Since we know the relation m 2 e . , then we can expect |A fl | \A L \. The A L,R 

contain the contribution of the neutralino loop and the chargino loop as seen in fig. 1 . 

We calculate the branching ratio using (|3.2|) and the formulas in Appendix C. In 
order to clarify parameter dependence, let us present an approximate estimation. The 
decay amplitude is approximated as 


\A R \ 2 ~ 


oi\ |(Am : 


L)iJ 1 


167T 2 


m 


tan 2 /3 , 


(3.3) 


where a 2 is the gauge coupling constant of SU(2) L and ms is a SUSY particle mass. 
The RGE’s develop the off-diagonal elements of the slepton mass matrix and A-term. 
These terms at the low energy are approximated as 


(Am 


LJ13 


(6 + 2a5 )^ (Y , YJj . ln V. 


167T 2 


m r 


(3.4) 


where M x is the GUT scale. Therefore, off-diagonal elements of (Y^Y u )ij are the 
crucial quantity to estimate the branching ratio. 
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As discussed in section 2, (Y^Y^- is given by neutrino masses and mixings at 
the electroweak scale. Therefore, we can compare the quantity (YjY„) zj among the 
cases of three neutrino mass spectra: the degenerate, the inverse-hierarchical and the 
hierarchical masses. In this section, we present numerical results in these three cases. 
Here, we use eq.( [3.2| ) and the vertex functions in Appendix C for the calculation of the 
branching ratio including the RGE’s effect. 


3.1 /i —» e + 7 


We present a qualitative discussion on (YtY„) 21 before predicting the branching ratio 
BR(/r —> e+y). This is given in terms of neutrino masses and mixings at the electroweak 
scale as follows: 

(YIY u ) 2 i = ^ [U^U^m^ - rriv i) + Uf, 3 U* 3 (m v3 - m vl )] , (3.5) 

W 

where v u = usin/3 with v = 174GeV is taken as an usual notation and the unitarity 
condition of the lepton mixing matrix elements is used. Taking the three cases of 
the neutrino mass spectra: the degenerate, the inverse-hierarchical and the normal 
hierarchical masses, one obtains the following forms, repectively, 


(Y$Y v ) 2 i - 


M R Am ltm 

V2v?. 2 m u 




y/2 e2 &ml 


e3 


atm 


(Degenerate) 


Mn 

V2v 


V Am atm 


1 tt* Am l 
U £>' 


2y/2 Am 2 


- U, 


e3 


atm 


(Inverse) 


(3.6) 


Mn 

V2v 


V Am atm 


--UI 


Am| + 


e3 


(Hierarchy) 


y/2 e2 ^Am 2 tm 

where we take the maximal mixing for the atmospheric neutrinos. Since U e2 — l/v^ 
for the bi-maximal mixing matrix, the first terms in the square brackets of the right 
hand sides of eqs.([h6|) can be estimated by putting the experimental data. For the 
case of the degenerate neutrino masses, (YjY^i depends on the unknown neutrino 
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mass scale m u . As one takes the smaller m„, one predicts the larger branching ratio. 
In our calculation, we take m v = 0.3eV [] , which is close to the upper bound from the 
neutrinoless double beta decay experiment |25[ , and also leads to the smallest branching 
ratio. 

We also note that the degenerate case gives the smallest branching ratio BR(/i —> 


e + 7) among the three cases as seen in eqs. (|3.6| ) owing to the scale of m„. It 
is easy to see the fact that the second terms in eqs.( |3.6| ) are dominant as far as 
I/ e3 ~ 0.01 (degenerate), O.Ol(inverse) and 0.07(hierarchy), respectively. The magnitude 
and the phase of I/ e3 are important in the comparison between cases of the inverse- 
hierarchical and the normal hierarchical masses. In the limit of I/ e3 = 0, the predicted 
branching ratio in the case of the normal hierarchical masses is larger than the other 
one. However, for C/ e3 — 0.2 the predicted branching ratios are almost the same in 
both cases. 

At first, we present numerical results in the case of the degenerate neutrino masses 
assuming M R = M R 1. The magnitude of M R is constained considerably if we impose 
the b — t unhcation of Yukawa couplings |26|]. In the case of tan (3 < 30, the lower 
bound of M r is approximately 10 12 GeV. We take also M R < 10 14 GeV, in order that 
neutrino Yukawa couplings remain below 0(1). Therefore, we use M R = 10 12 , 10 14 GeV 
in our following calculation. 

We take a universal scalar mass (mo) for all scalars and ao = 0 as a universal 
A-term at the GUT scale (M x = 2 x 10 16 GeV). The branching ratio of /r —> e + 7 
is given versus the left-handed selectron mass for each tan/3 = 3, 10, 30 and 
a fixed wino mass M 2 at the electroweak scale. In hg.^, the branching ratios are 
shown for M 2 = 150, 300 GeV in the case of I/ e3 = 0.2 with Mr = 10 14 GeV and 
m u = 0.3eV, in which the solid curves correspond to M 2 = 150GeV and the dashed 


1 Recently, a positive observation of the neutrinoless double beta decay was reported in where 
the degenerate neutrino mass of m v = 0.3eV is a typical one. 


11 












Figure 2: Predicted branching ratio BR(p —> e + 7 ) versus the left-handed selectron 
mass for tan/? = 3, 10, 30 in the case of the degenerate neutrino masses. Here M R = 
10 14 GeV and f / e 3 = 0.2 are taken. The solid curves correspond to M 2 = 150GeV and 
the dashed ones to M 2 = 300GeV. A horizontal dotted line denotes the experimental 
upper bound. 



Figure 3: Predicted branching ratio BR(p —> e + 7 ) versus the left-handed selectron 
mass for tan /3 = 3, 10, 30 in the case of the degenerate neutrino masses. Here 
M r = 10 12 GeV and U e 3 = 0.2 are taken. The solid curves correspond to M 2 = 150GeV 
and the dashed ones to M 2 = 300GeV. 
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Inverse-hierarchical type MR=10 14 GeV 



Figure 4: Predicted branching ratio BR(p —> e + 7 ) versus the left-handed selectron 
mass for tan/? = 3, 10, 30 in the case of the inverse-hierarchical neutrino masses, ffere 
M r = 10 14 GeV and U e 3 = 0.2 are taken. The solid curves correspond to M 2 = 150GeV 
and the dashed ones to M 2 = 300GeV. 


ones to M 2 = 300GeV. The threshold of the selectron mass is determined by the 
recent LEP 2 data J27J for M 2 = 150GeV, however, for M 2 = 300GeV, determined 
by the constraint that the left-handed slepton should be heavier than the neutralinos. 
As the tan (3 increases, the branching ratio increases because the decay amplitude 
from the SUSY diagrams is approximately proportional to tan/3 ]7|. It is found that 
the branching ratio is almost larger than the experimental upper bound in the case 
of M 2 = 150GeV. On the other hand, the predicted values are smaller than the 
experimental bound except for tan (3 = 30 in the case of M 2 = 300GeV. 

Our predictions depend on M R strongly, because the magnitude of the neutrino 
Yukawa coupling is determined by M R as seen in eq. ( P-5|) . If M R reduces to 10 12 GeV, 
the branching ratio becomes 10 4 times smaller since it is proportional to M\. The 
numerical result is shown in hg.|3]. We will examine a model [^, which gives the 
degenerate neutrino masses with U e 3 ~ 0.05 in section 4. 
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Inverse-hierarchical type MR = 10 14 GeV , U e 3=0.05 



Figure 5: Predicted branching ratio BR(p —> e + 7 ) versus the left-handed selectron 
mass for tan/? = 3, 10, 30 in the case of the inverse-hierarchical neutrino masses, ffere 
M r = 10 14 GeV and U e 3 = 0.05 are taken. The solid curves correspond to M 2 = 150GeV 
and the dashed ones to M 2 = 300GeV. 


Next we show results in the case of the inverse-hierarchical neutrino masses. As 
expected in eq. (|3.6|) , the branching ratio is much larger than the one in the degenerate 
case. In hg.||, the branching ratio is shown for M 2 = 150, 300 GeV in the case of 
U e3 = 0.2 with M r = 10 14 GeV. In fig. I, the branching ratio is shown for U e 3 = 0.05 
with M r = 10 14 GeV. The M R dependence is the same as the case of the quasi¬ 
degenerate neutrino masses. The predictions almost exceed the experimental bound 
as far as U e 3 > 0.05, tan/3 > 10 and M R ~ 10 14 GeV. This result is based on the 
assumption M R = M R 1, however, it is not guaranteed in the case of the inverse- 
hierarchical neutrino masses. We will examine a typical model [3(J, which gives M R 7 ^ 
M r 1 in section 4. 

For comparison, we show the branching ratio in the case of the hierarchical neutrino 
masses in fig.[]. It is similar to the case of the inverse-hierarchical neutrino masses. The 
branching ratio in the case of the degenerate neutrino masses is 10 2 times smaller than 
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Hierarchical type M R =10 14 GeV 



Figure 6 : Predicted branching ratio BR( / u —> e + 7 ) versus the left-handed selectron 
mass for tan/3 = 3, 10, 30 in the case of the hierarchical neutrino masses. Here 
M r = 10 14 GeV and U e3 = 0.2 are taken. The solid curves correspond to M 2 = 150GeV 
and the dashed ones to M 2 = 300GeV. 

the one in the inverse-hierarchical and the hierarchical neutrino spectra. 

In our numerical analyses we assumed ao = 0 at the GUT scale M x for simplicity. 
Let us comment on the A-terrn dependence, namely ao 7 ^ 0 at M x . We estimate 
the branching ratio for ao = ±1 at M x (A = Yaomo). In the degenerate type, the 
predicted branching ratio is 1.02(ao = 1), 1.07(ao = —1) times as large as the one in the 
case of a 0 = 0 (tan (3 = 30, U e3 = 0.2). In the inverse-hierarchical type, the predicted 
branching ratios are 1.56(ao = 1), 1.54(ao = —1) times as large as the one in the case 
of ao = 0 (tan/3 = 30, U e 3 = 0.2). Therefore the A-term dependence is insignificant in 
our analyses. 

In our calculations, we use the universality condition at M x . We also examine the 
no-scale condition mo = 0 at M x . It is found that the predicted branching ratio is 10 
times smaller than the one in the case of non-zero universal scalar mass. 
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3.2 Non-degeneracy Effect of 

The analyses in the previous section depend on the assumption of M Rl = M r2 = M R3 = 
Mr. In the case of the quasi-degenarate neutrino masses in eq. (|2.12|) this complete 
degeneracy of may be deviated in the following magnitude without fine-tuning: 


M\ 


r3 


m 2 r1 


~ 1 ± 


Am 


atm 


Mi 


R2 


mf 


M 2 


~ 1 ± 


Ami 


o 


Rl 


mf. 


(3.7) 


Therefore, we parametrize M R as 


= M r 


( 1 0 0 

0 1 + e 2 0 


(3,8) 


0 0 1 + £3 

where e 2 ~ Am 2 & /2m 2 and e 3 ~ Am 2 tm /2ml. By using eq. (|2.6|) , we obtain 


Y t Y _ Majj 

L 1/ L V — o ^MNS 

v z 

'-'ll 


/ y/rfhil 

0 

0 N 

\ 

( y/m^l 

0 

0 

0 

V m v 2 

0 

K 

0 

y/m v2 

0 

V 0 

0 

a/^3 , 

/ 

l 0 

0 

) 


where 


K = R f 


/ 1 0 0 

0 1 T e 2 0 I R 

V 0 0 1 + e 3 


Then, we have 


M 3 

(Y^Y„) 2 i = , 

Vu i,j 


with 


U 


MNS 1 


(3.9) 


(3.10) 


(3.11) 


Kij — Sij + E 2 R 2 i R 2j + £ 3 R:>,,R:>, : j , 
where we used RdR = 10. So, we get 


YtY „) 21 = (YlY; 


V 21 


Mr.0, 


A(Y*Y„) 21 , 


2 we assume R to be real for simplicity. 


(3.12) 


(3.13) 
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where the first term is the (YtY v ) 21 element in eq.( |3.5|) , which corresponds to the 
M roc1 , while the second term stands for the deviation from it as follows: 


^O^Y„) 2 i — —y ^ U 2l U 1 j \/fflm^Ul v j (e 2 R' 2 iR' 2 j + £3^31^3^ ■ (3-14) 

V u i,j 

In order to estimate the second term, we use £ 2 = 0.0001 and £3 = 0.01 taking account 
of e 2 — Am 2 /2m 2 and £3 ~ Am^ tm /2mJ, where m u = 0.3eV is put. Since rn vl ~ m V j 
and Rij < 1 , we get 

M 3 

A(Y,tY „) 21 ~ ~?Z U 2 iUi j 7 n u £ 3 R 3 i R 3j 
Vu i,j 

~ ^ 27 2 77lu£ 3 ~ 3 ' 5 X 10 ~ 3 ' (3 ' 15) 


Taking this maximal value, we can estimate the branching ratio as follows: 

5i?(non-degenerate M R ) < /2.6 + 3.5 \ 2 
Hi?(degenerate M R ) ~ V 2.6 / 


Therefore, the enhancement due to the second term is at most factor 5. This conclusion 
does not depend on the specific form of R 

Consider the case of the inverse-hierarchical type of neutrino masses. We take 
£ 2 ~ 0.01 with the similar argument of the quasi-degenerate type neutrino masses, 
because m u \ and m v2 are almost degenerate and £ 2 — Am 2 / 2 Am 2 tm in this case. 
Then, we get 


A(Y*YJ 21 = U 2 iUijm u i£ 2 R2iR2j 

V u i,j 

M r 1 _ 9 

< —-- -=m v i£ 2 ~ 0.063 x 10 " , 

K 2V2 

where we assume £ 2 > £3 and use m v3 ~ 0, m v 1 ~ m v2 — 0.054eV and R t] < 1. 
the maximal value, we get 


(3.17) 

Taking 


Ri?(non-degenerate M R ) /2.7 + 0.063\ 2 

-< - ~ 1.04 

5i?(degenerate M R ) — V 2.7 J 


(3.18) 
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Degenerate type 


M R = 10 14 GeV , U e3 =0.2 


1CT 6 

nr 7 

1CT 8 

BR(t -> py) nr 9 
1CT 10 
1CT 11 
1CT 12 

1(T 13 1CT 12 1CT 11 1CT 10 1CT 9 1CT 8 1CT 7 1CT 6 
BR(u^ey) 

Figure 7: Predicted branching ratio BR(t —> /i + 7 ) versus BR(/r —> e + 7 ) for tan (3 = 
3, 10, 30 in the case of the degenerate neutrino masses. Here M R = 10 14 GeV, m v = 
0.3eV and U e3 = 0.2 are taken and the left-handed selectron mass is taken as same as 
in figs.2-4. 

Thus, the effect of the A(YjY „) 21 is very small in the case of the inverse hierarchical 
neutrino masses. 

These discussions in this subsection are also available qualitatively for the r —> /i +7 
process. 



3.3 r —» fj, + 7 

Let us study the r —> /1 + 7 process. In this case, we should discuss 


M, 


(YjY „) 32 = —Y UnU^m^ - m vl ) + U T3 U* 3 (m„ 3 - m„i) 


(3.19) 


It should be stressed that it is independent of U e3 in contrast to (yJy„) 21 . Therefore 
we can determine the following form of (YjY ^) 32 at the electroweak scale by using the 
bi-maximal mixing matrix: 
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Inverse-hierarchical type 


M R =10 14 GeV , U e 3=0.2 


1 CT 6 

10- 7 

1 CT 8 

BR(x -» py) 1CT 9 
1 (T 10 
1 CT 11 
1 CT 12 



10 13 10 12 10" 11 icr 10 1CT 9 1CT 8 icr 7 1CT 6 


BR(p->ey) 


Figure 8 : Predicted branching ratio BR(r —> fi + 7 ) versus BR(/i —> e + 7 ) for tan f3 = 
3, 10, 30 in the case of the inverse-hierarchical neutrino masses. Here M R = 10 14 GeV 
and U e 3 = 0.2 are taken and the left-handed selectron mass is taken as same as in 
figs. 2-4. 


YtY „) 32 ~ ^ 


l Am 2 & + i Am atm 


4 2 m h 


2 2 m v 


M„ Ami 


atm 


Avl m v 


, (Degenerate) 


M, 


1 A ”i -iJAml 


8 \/ A "4 


atm 


atm 


Mn 
’2 vl 


\/AmL , (Inverse) 


(3.20) 


Mj 

vl 


■|V Am 0+ 2 V Am atm 


~ JAm 2 

- 2v 2 V 


. (Hierarchy) 


We see that the case of the inverse-hierarchical masses and the hierarchical masses are 
almost the same as seen in eqs.( |3.20| ). 

Let us present numerical results of BR(r —> // + 7 ) |HJ versus BR(/i —> e + 7 ) JTI 


in the case of the degenerate neutrino mass, in which tan/3 = 3, 10, 30 are taken. 
In hg.[7j. the branching ratio is plotted for M 2 = 150, 300 GeV for I / e3 = 0.2 with 
M r = 10 14 GeV. Dotted lines are the experimantal upper bounds for BR(r —> /x + 7 ) 
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and BR(/i —> e + 7 ), respectively. The dependence of tan/3 is the same as the case 
of ft —> e + 7 . It is found that the branching ratio is completely smaller than the 
experimental upper bound in the case of r —> p + 7 in contrast with the case of 
p —> e + 7 . 

Next we show the results in the case of the inverse-hierarchical neutrino masses. 
As expected in eqs.( |3.20| ), the branching ratio is much larger than the one in the 
degenerate case. In fig.[^, the branching ratio is shown for ilf 2 = 150, 300 GeV in the 
case of U e 3 = 0.2 with M R = 10 14 GeV. In conclusion, the predicted branching ratio is 
larger than the one in the case of the degenerate neutrino mass, and it is almost smaller 
than the experimental upper bound for r —*■ p + 7 in contrast with /1 —> e + 7 . The 
constraint of BR(p —> e+ 7 ) is always severer than the one in the case of BR(r —> /i+ 7 ). 


4 Typical models and numerical analyses 

4.1 Ss L x S% R flavor symmetry model - Degenarate type 


In this section we examine the neutrino model proposed by Fukugida, Tanimoto and 
Yanagida |28|, which derives the quasi-degenerate masses, rn v \ ~ m V 2 ~ m„. 3 . This 
model is based on the S$ L x Ss R flavor symmetry |32|. Taking M R = M R 1, the neutrino 
Yukawa coupling is given as follows: 


(4.1) 




( 1 0 0 \ 

O 

O 

O 

Y, = Y u0 


010 + 

0 e u 0 



0 1/ 

^ 0 0 6„ ) 


where we take the diagonal basis for the neutrino sector. The first matrix is the S;$ L in¬ 
variant one, and the second one is the symmetry breaking term. The parameters Y u 0 , 
and S u are constrained by the experimental values of Am^ tm and A m 2 Q . Therefore, the 
flavor mixings come from the charged lepton Yukawa couplings. 

The charged lepton Yukawa coupling is given by the symmetry breaking parameters 
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Figure 9: Predicted branching ratio BR(p —> e + 7 ) versus the left-handed selectron 
mass for tan/? = 3, 10, 30 in the case of the S 3 l x S 3 r flavor symmetry model, ffere 
M r = 10 14 GeV are taken. The solid curves correspond to M 2 = 150GeV and the 
dashed ones to M 2 = 300GeV. 

€ 1 , 5i as follows: 




f 1 

1 

1 N 


( ~ ei 

0 

° ^1 

Y e = Y e0 


1 

1 

1 

+ 

0 

+Q 

° 



l 1 

1 

1? 

/ 

l 0 

0 

+5/ )\ 


(4.2) 


Since Y e0 , e; and 5i are fixed by the charged lepton masses, one gets the lepton mixing 
matrix elements as follows : 


U 


MNS 




/ l/y/2 -l/y/2 ^2/3^/mJm^ \ 

1/V6 1/Vs -2/V6 
V 1/V3 1A/3 l/y/3 ] 


(4.3) 


As a result, we see U e3 = ^2/Z\Jm e /m ll ~ 0.05 from eq. (|4.3|) . 

We estimated the branching ratio of the processes p —> e + 7 and r —> p + 7 by 
using U e3 = 0.05. We show the branching ratio for M 2 = 150 GeV and 300GeV taking 
tan/? = 3,10, 30 in fig|]. Because of the smallness of U e3 , we see that BR(p —> e + 7 ) is 
smaller than the experimental upper bound except for tan/? = 30, and M 2 = 150GeV. 
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We have also estimated the branching ratio BR(r —> p + 7 ) for tan/3 = 30, which 
is much smaller than the experimental bound BR(t —> fi + 7 ) < 1.1 x 10~ 6 . Thus, the 
p ^ e + 7 process provides the severe constraint compared with the r —> fi + 7 in the 
present experimental situation. 


4.2 The Shafi-Tavartkiladze model - Inverse-hierarchical type 


The typical model of the inverse-hierarchical neutrino masses is the Zee model [p3j| , 
in which the right-handed neutrinos do not exist. However, one can also consider a 
Yukawa texture which leads to the inverse-hierarchical masses through the see-saw 
mechanism, namely the Shafi-Tavartkiladze model [^[ . 

Shah and Tavartkiladze utilize the anomalous U (1) flavor symmetry [33]. In this 
model, due to the Froggatt-Nielsen mechanism |F|, one of the Yukawa interaction term 
in the effective theory is given as 

/ qi \ m ij 


e c R *LjH d 


M, 


Pb 


(4.4) 


where e c Ri and Lj are the right-handed charged lepton and the left-handed lepton dou¬ 
blet, respectively, H d is Higgs doublet, and S is singlet field. The effective Yukawa 
couplings are given in terms of 


A = ~ 0.2 


Mr 


(4.5) 


pi 


The neutrino mass matrix is given in Appendix D. Fixing the [/(1) flavor charges 
/c, n, k' as k= 0, ?r=2, k'= 2, which is consistent with neutrino mass data, the Yukawa 
coupling is given as 


Y„ = 


A 4 A 2 A 2 
1 0 0 


(4.6) 


and the right-handed neutrino Majorana mass matrix is given as 

M r = Mr 


( A 4 1 
lR * 1 0 


22 


(4.7) 













Shafi-Tavartkiladze model 



Figure 10: Predicted branching ratio BR(/i —» e + 7 ) versus the left-handed selectron 
mass for tan/3 = 3, 10, 30 in the case of the Shafi-Tavartkiladze model[j30fl. The solid 
curves correspond to M 2 = 150GeV and the dashed ones to M 2 = 300GeV. 


It is remarked that right-handed neutrinos contain only two generations in this model. 
In eq. (|4.6|) , components 2-2 and 2-3 must be zero for the sake of holomorphy of super¬ 
potential, it is called SUSY zero. The neutrino mass matrix is given by the see-saw 
mechanism as 


m„ = YjM-'Y^ = 


2„,2 / A 2 1 1 


\\ 

Mr 


1 0 0 
V 1 0 0 


(4.8) 


where the order one coefficient in front of each entry is neglected. This mass matrix 
gives the inverse-hierarchical neutrino masses. The Yj Y u is given as 


/ 1 + A 8 A 6 A 6 


Y+Y — 

x V x ^ 


A 6 

V A« 


A 4 A 4 
A 4 A 4 


It is noticed that the component (YjY 1/ ) 2 i is suppressed as 


(4.9) 


(YiY„) 2 i ~ A 6 ~ 0( IQ' 5 ) 


(4.10) 
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Then, we expect that the branching ratio of /j —> e + 7 in this model is much smaller 
than the one in the case of (M R ) 3x3 = M R ( l ) 3x3 in section 3. 

In fig.|lO], the branching ratio is shown for M 2 = 150, 300 GeV. The predictions 
are given by taking A = 0.2 and all of order one coefficients in the Yukawa couplings 
are fixed to be one. The predicted value is much smaller than the one in the inverse- 
hierarchical case discussed in the section 3. Because (Y^Y^i is proportional to A 6 , 
the smallness of the branching ratio is understandable. 

5 Summary and Discussions 

We have investigated the lepton flavor violating processes p —> e + 7 and r —> /1 + 7 , 
in the framework of the MSSM with the right-handed neutrinos. Even if we impose 
the universality condition for the soft scalar masses and A-terms at the GUT scale, 
off-diagonal elements of the left-handed slepton mass matrix are generated through 
the RGE’s running effects from the GUT scale to the right-handed neutrino mass scale 
M r . We have taken the LMA-MSW solution for the neutrino masses and mixings. 

The branching ratios of /1 —> e + 7 and r —> p + 7 processes are proportional 
to KYfYy).J 2 . Since (Y^Yj,)^ depends on the mass spectrum of neutrinos, we can 
compare the branching ratio of three cases of neutrino mass spectra: the degenerate, 
the inverse-hierarchical and the hierarchical case. 

First, we have studied the three types in the case of p —> e + 7 , in which we 
take M r = M It l. For the case of the degenerate neutrino masses, the branching 
ratio depends on the unknown neutrino mass m u . We have taken m„ = 0.3eV, which 
gives us the largest branching ratio. It is emphasized that the magnitude of U e3 is 
one of important ingredients to predict BR(p —> e + 7 ). The branching ratio of the 
inverse-hierarchical case almost exceeds the experimental upper bound and is much 
larger than the degenerate case for M 2 = 150GeV and M 2 = 300GeV. In general, we 
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expect the relation BR,(degenerate) <CBR, (inverse-hierarchical) < BR(hierarchical). The 
effect of the deviation from M R = M R 1 has been estimated. The enhancement of the 
branching ratios are at most factor five in the case of the quasi-degenerate neutrino 
mass spectrum. 

Second, we have studied the three cases in r —> /x + 7. It is noticed that branching 
ratio is independent of C / e3 in contrast to the case of n —> e + 7. For the degenerate 
neutrino masses, the branching ratio is completely smaller than the experimental upper 
bound. For the inverse-hierarchical neutrino masses, the branching ratio is almost 
smaller than the experimental bound. The constraint of BR(yU —» e + 7) is always 
severer than the one in the case of BR(r —> fj, + 7 ). 

Finally, we have investigated the branching ratio of fj, —> e + 7 in the typical models 
of the degenetate and the inverse-hierarchical cases. Since the 63 L x ,S' 3/i model, which 
is a typical one of the degenerate case, predicts U e 3 ~ 0.05, the branching ratio is much 
smaller than the case of 7 / e3 ~ 0.2. The Shafi-Tavartkiladze model, which is a typical 
one of the inverse-hierarchical case, predicts the very small branching ratio. Thus, the 
models can be tested by the fi —> e + 7 process. 

The branching ratio of fi —> e + 7 and r —> /x + 7 will be improved to the lebel 10 “ 14 
in the PSI and 10 ~b~ 8 ) j n the B factories in KEK and SLAC, respectively. Therefore, 
future experiments can probe the framework for the neutrino masses. 
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A Yukawa Matrix 


The Yukawa matrix is determined in general as follows 
mass matrix is given as 


101. The left-handed neutrino 


m u = (Y u v u ) T M r 1 (Y u v u ) , 


(A.l) 


via the see-saw mechanism, where v u is the vacuum expectation value (VEV) of Higgs 
H u . One can always take the diagonal form of the right-handed Majorana neutrino 
mass matrix M R = M R ag . The neutrino mass matrix m v is diagonalized by a single 
unitary matrix 

m v ias = U^hi.Umns , (A.2) 


where U MNS is the MNS matrix, 
square roots 


In eqs.( A.l ) and ( A.2 ), one can divide M R ag into 


m 


diag 

V 


ULsYj (M“T‘ Y„U„«f2 
Um NS YJv^MST' sj . 


(A.3) 


Multiplying the inverse square root of the matrix m)( iag from both right and left hand 
sides of eq.( |A-3|) , one gets the following form 


m diag')- 1 pj T Y t 

^ MNS A V 


(M R as ) \? u 




rn 


diag 


,-l 


- dT 


R R , 


(A.4) 


where one has defined the following complex orthogonal 3x3 matrix 


R = v u \J (M R as ) Y^UmnsV (mf* 


,-i 


(A.5) 


and R depends on models. Therefore, one can write the neutrino Yukawa coupling as 

1 C3Z. n —... (A.6) 


Y„ = — x/MS- R ,/m#- U 


V MNS ’ 
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or explicitly 



0 

o ^ 


/ 

0 

° ^ T 


0 

V Mr2 

0 

R 

0 

y/m v2 

o U MNS . 

(A.7) 

V o 

0 

V Mr 3 ) 


l o 

0 

■sjniu-s ) 



B RGEs 

B.l From M x to M R 


d/i 

n—Y ij 
^d/J S 


n—Y ij 

l 1 i *- V 

a/i 


n—Y l i 
^dn u 


n—Y ij 
^dfi d 


1 

—b ig l, (6 1 ,6 2 ,6 3 ) = (y, 1,-3), 

b■ (7? 

= -^ (* = 1, 2, 3) , 

27T 47T 


9 


0? - 3(^2 + 3Tr(Y d Yl) + Tr(Y e Yt) 1 Y? 


167T 2 LI 5 
+3(Y e Y]Y e )^' + (Y e YlY u ) ij 


-^gl ~ 3s| + 3Tr(Y u Y,t) + Tr(Y.Yt) [> Yj? 


167T 2 

+3(Y„YjY„)« + (Y„Y*Y e )« 


^9? - 39? - f 9f + 3IY(Y„Y*) + Tr(Y„Yj) Yj 


167T 2 

+3(Y (t Y^Y. u )* J + (Y tt YjY d )« 

7 16 

-Jfi ~ 3gl - ys! + 3Tr(Y„Y2) + Tr(Y e Y|) \ Y‘j 


167T 2 


+3(Y d YjY d ; 


y 


(Y d YtY tt )«" 


/i 


A 

d/j, 


m 


2\j 


a/i J 




167T 2 


(m 2 £ Y|Y c + YtY.m?) 2 + (m?Y»Y„ + Y*Y„m 


It 


+2 (Y*m?Y e + my Y|Y e + A*A.) 
+2 (Yjm|Y„ + my Y*Y„ + A 2 A 

- (|s? Mi 1 2 + e 9 f |m 2 | 2 ) y], 

2 (m 2 Y e Yj + Y.Yjmi) 1 


167T 2 


+4 (Y e m?Y* + my Y e Y* + A,A*)’ - - 9 \\M,\ 2 S) 
2(m 2 Y„Yt+Y„Ytm 2 )’ 


I6n 2 


32 



+4 (Y,m? Yt + m 2 Hu Y u Yl + A„At)‘] , 

= I ^2 [{-^i - + 3Tr(YlY d ) + Tr(YtY e )} A« 

+2 {-plM, - 3g 2 2 M 2 + 3Tr(YlA d ) + Tr(Y|A e )} Y” 
+4(Y e YtA e )« + 5(A c YtY e )tf + 2(Y e Yt A v f + (AgY^Y^)*- 7 ] , 

= 1^2 3 ^2 + 3Tr ( Y i Y «) + Tr ( Y i Y ^)} A t J 

+2 - 3g 2 M 2 + 3Tr(Yt A tt ) + Tr(Y^A v )| Y? 

+4(Y I/ YtA„)« + 5(A,YtY,)^' + 2(Y I/ Y+A e )« + (A,Y]Y e )^] , 
4 A “ = I^[{-l^- 3 ^-fft 2 + 3TV( Y .Y„) + IV(YlY„)}A« 

+2 {-4|s?Afi - 3 s|M 2 - ^sfM 3 + 3Ti'(Y*A„) + Tr(Y,U„)} Y« 
+4(Y„Y*A„)« + 5(A„YiY„)« + 2(Y,Y*Aj)« + (A„Y»Yj)«] , 

= I^[{4 9 ?- 392 -f 92 + 3Tr(Y ^ ) + TV(Y Y)}A« 

+2 - 3» 2 2 M 2 - YsfM 3 + 3Tr(Y'A,,) + Tr(Y|A e ) j Y« 

+4(Y J Y»A tl )« + 5(A,,Y»Y,,) W + 2(Y d Y*A„)« + (A,Y*Y„)«] , 

d -J^ 

/^ (^ = 16^2 [ 6Tr ( m l Y " Y « + Y “( m “ + )' Y « + A l A «) 

+2Tr (mfY+Y,, + Y^(m? + m 2 Hu )Y v + A^A,,) 

— ^5 , i|-^i| 2 + 6(^21A^21, 

^ [ 6Tl ' ( m | Y rf Y d + Y l( m d + m H d ) Y d + A 1 A rf) 

+2Tr (m?YtY e + Yj(m \ + ro^) Y e + A+A*) 

— + 6^?2 I AY 2 | 2 ^) • 
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C.l Mass matrix and mixings 


In this appendix, we give our notation of SUSY particle masses and mixings in our 
calculation. 

The slepton mass M 2 term is 


with 


(C4) 


m, 

m 2 



( m Dp 


i m 




[ III 


lrJv 


(m l)ij + m 2 e .8ij + m 2 z 8ij cos 2/3(-- + sin 2 9 W ) , 

(m e R )ij + m 2 ei 8ij - m\8ij cos 2/3 sin 2 8 W , 

A.®- v cos/ 3 


V2 


m ei ntan/3 , 


(C.l) 


(C.2) 

(C.3) 

(C.4) 


where (m 2 )y and (m 2 ) y are 3x3 matrix. The slepton mass matrix can be diagonalized 


as 


I/lM 2 !/! 1 = (diagonal) , 


(C.5) 


where U* is a real orthogonal 6x6 matrix. 
The chargino mass term is 


£ = (wr^H^) 


M 2 \/2m w cos /3 
y/2m w sin /3 // 


H 


+ h.c. 


(C.6) 


1 L 


The chargino mass matrix can be diagonalized as 


O r M c OJ = diag(M--, M-~) , 


(C.7) 


where O l and O r are real orthogonal 2x2 matrix. The mass eigenstate Xal and \ AR 
{A = 1, 2) are 


XlL 
X.2 L 


= 0 L 




f XlR 

\ X.2 R 


= o f 


h 2r 



(C.8) 
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and 


Xa Xal Xar 


(-4=1,2) 


(C.9) 


forms Dirac fermion with mass M-~. 

Xa 

The neutralino mass term is 


t B l \ 


-C = 


-(B l ,W° l ,H 0 1l ,H° 2l )M n 


W L U 

Hj L 

\ h° 2l J 


+ h.c. , 


(C.10) 


where 


M n = 


( M 1 0 

0 m 2 

—rriz sin d w cos (3 mz cos cos /3 

V mz sin 0 W sin (3 —mz cos 6yj sin (3 


—rriz sin $w cos (3 mz sin $w sin (3 ^ 
mz cos cos /3 —mz cos 6 w sin /3 

0 — n 

-h 0 / 


(C.ll) 


The neutralino mass matrix can be diagonalized as 


O n M n O t = (diagonal) , 


(C.12) 


where O n is a real 4x4 orthogonal matrix. The mass eigenstate is 


Xal — {Os) 


abXbl = l,---,4:) 


X° BL = (B h , W, 


o uO 


TT U TT U \ 

11 2l) 


(C.13) 


and 


Xa — Xal + Xar — 1, * * *, 4) 


forms Mojorana spinor with mass M^o. 
The chargino vertex functions are 


^ eAX 

^ eAX 


— 92(0r)aiUx,\ 4 


= 9 2 


rn f 


y/2m w cos /3 


(OlUuz* , 


(C.14) 


(C.15) 

(C.16) 


35 



and the neutralino vertex functions are 


= --%([— (0 N ) A2 -(0 N ) Al tane yr \U t Xil + me J O n ) a 3 U x A , (C.17) 
V 2 I cos p j 


9 2 




= 

v/2 1 m w cos/3 


(On) A'iU x a 2(O n ) a i tan 9 w U x 4 


(C.18) 


C.2 Decay amplitudes A 


L.R 


For the amplitudes A L ’ R , there are the contributions of the chargino loop and the 
neutralino loop: 


J^L,R _ J^(c)L,R _|_ n)L,R 


(C.19) 


The contributions from the chagino loop are 




vx 


(~i l(£) 


32n 2 m 2 „ | JAX lAX 6(1 — x AX ) 4 


^R(e)* ^x A 1 

' ^jAX^iAX /-i \q 

nrij (1 - x Ax y 


(2 + 3x AX - 6x 2 AX + x AX + 6x AX In x AX ) 
(-3 + 4x AX - x 2 AX - 2 lnx AX ) , (C.20) 


jy°) R — a(°) l I 


(C.21) 


where x AX is defined as 


ML 
_ x A 

%AX — 9 

m lx 


(C.22) 


Here rn ^„ is the sneutrino mass and M-- is the chargino mass. The contributions from 


the neutralino loop are 


lx 


A{U)L = 32L mL” N iAx N i!ix * _ yA j4 (1 ~ 6 Vax + 3 y\ x + 2 y\ x - 6 y 2 AX In y AX ) 


_l_ ]\f R (l)* _jj_ 

' I ' l jAX ly iAX \q ' 

m (1 - y Ax y 


(1 + Vax + AX In y AX ) 


A ( n )« _ a^ l 


L<—>R i 


, (C.23) 
(C.24) 


where y AX is defined as 


Vax — 


Ml o 

Xa 


m 


(C.25) 




Here is the charged slepton mass and M^o is the neutralino mass. 
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D Anomalous U{ 1) flavor symmetry 


We review the anomalous [/(1) flavor symmetry |I4| which is utilized in the Shafi- 
Tavartkiladze model |3U] discussed in section 4. The anomalous 17(1) flavor symmetry 
can arise from string theory. The cancellation of this anomaly is due to the Green- 
Schwarz mechanism |36|. The associated Fayet-Iliopoulos term is given as [37 


£ j d 4 9V A with 


£ = 


lA 1 V 1 p\ 

192vr 2 


TrQ . 


(D.l) 


The D-term is given as 


t d ’ = t (EQ.Irf + e) 2 . 


9 2 a 


(D.2) 


where Q a is the ‘anomalous’ charge of ip a . For 17(1) breaking, we introduce the singlet 
field S under the SM gauge group with 17(1) charge Q s . Assuming Tr Q > 0, we can 
ensure the cancellation of D A in ( |D.2j ). Taking Q s = —1, we can ensure the non-zero 
VEV of S: ( S ) which is given as ( S) = ^. 

Due to the Froggatt-Nielsen mechanism, Yukawa interaction term in the effective 
theory is given as 


e c RJ L t H d 



(D.3) 


where e c R j and Li are the right-handed charged lepton and left-handed lepton doublet, 
respectively, H d is Higgs doublet, and S is singlet field. The effective Yukawa couplings 
are given in terms of 


M p i 


(D.4) 


In order to make the interaction term eq. (p.3| ) neutral, Shaft and Tavartkiladze 
assigned 17(1) flavor charges as follows: 


Qli — k + n , Q L2 — Q L3 — k , Q Nl — — Qn 2 — k + k , 
Qh u = QH d =0, Q s = -1 , 


(D.5) 
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where k,n, k' >0, n>k'. Then, they obtained 


Y„ = 


y2k-\-n-\-k y2 k+k y2k+k 


A 


n—k 


o 


o 


and 


Mr = M r 


y2k+2k y 

1 o 


In conculsion, the neutrino mass matrix is given as 

/ A n 1 1 


m = Y T M _1 Y v 2 = 

1 V 1VJ -R V U U 


\ 2k+n v 2 


1 0 0 
V i o o 


(D.6) 


(D.7) 


(D.8) 
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